Preprint TTK- 10-23 



Boltzmann hierarchy for the cosmic microwave background at second order 

including photon polarization 

M. Beneke and C. Fidler 

Institut fiir Theoretische Teilchenphysik und Kosmologie, 
RWTH Aachen University, D - 52056 Aachen, Germany 
(Dated: March 10, 2010) 

^— ^ ' Non-gaussianity and _B-mode polarization are particularly interesting features of the cosmic mi- 

, crowave background, as - at least in the standard model of cosmology - their only sources to first 

' order in cosmological perturbation theory are primordial, possibly generated during inflation. If the 

^ , primordial sources are small, the question arises how large is the non-gaussianity and B-mode back- 

■ ground induced in second-order from the initially gaussian and scalar perturbations. In this paper 

' we derive the Boltzmann hierarchy for the microwave background photon phase-space distributions 

. at second order in cosmological perturbation theory including the complete polarization informa- 

^\ ' tion, providing the basis for further numerical studies. As an aside we note that the second-order 

collision term contains new sources of B-mode polarization and that no polarization persists in the 
' tight-coupling limit. 

^ ! PACS numbers: 98.80.-k 

^ : 

I ■ I. INTRODUCTION 

o ■ 

ts: . . . . 

j/j ' The anisotropics of the cosmic microwave background provide an abundant source of information on the early 
history of the universe unrivaled in precision. Since most of the photons originate from the time of decoupling, 
I when the inhomogeneities of the universe were small, the anisotropies should be well described in linear perturba- 
^-H ■ tion theory around the Friedmann- Robertson- Walker background. Temperature anisotropies [l|, 0| and E-mode 
^ , polarization anisotropies [l|, H, |1] have been detected, and are found in agreement with the Standard Model, 
in which the source of the anisotropies consists of a gaussian, adiabatic and nearly scale-invariant spectrum of 
. primordial density perturbations. 
' The polarization pattern of the background radiation is of great interest for the following reason. In contrast 

fT^ ■ to the E mode, _B-mode polarization is not sourced by scalar density perturbations in the linear order. Thus, a 
I detection of i?-mode polarization would point directly to primordial vector or, more likely, tensor fluctuations 
' (gravitational waves) in a very early phase of the cosmological evolution. So far, however, i?-mode polarization 
• • , has not been observed, which together with the shape of the temperature perturbation spectrum indicates some 
. ^ ' suppression of tensor relative to scalar perturbations. Similarly, deviations from gaussian statistics constrain 
^ . inflation models and are therefore intensively investigated (see, e.g., the reviews 0,01). 

^ I The absence of non-gaussianity and B-mode polarization when the primordial fluctuations are purely gaussian 

- - ' and scalar holds, however, only in linear perturbation theory. Thus, if a small non-gaussian or i3-mode signal 
is observed, the question arises whether its origin is truly primordial, or whether it might be a second-order 
effect. While such an effect would naturally be expected at tensor-to-scalar ratios of order 10^^ (or /nl ~ 1 
for non-gaussianity), which is the size of perturbations in the microwave background, only a full second-order 
calculation can tell whether there are no enhancements. Such enhancements can reach a level relevant to 
observations, since the planned CMBPol experiment is sensitive to tensor-to-scalar ratios of order lO^'^ 
Indeed, several second-order sources of i3-mode polarization are already known. The most important is the 
weak-lensing effect, reviewed in 0], which converts E polarization to B polarization as the photons travel 
through the inhomogeneous universe 0] . The inhomogeneities and the i?-mode polarization are both at least of 
first order, so the resulting effect is at least of second order in perturbation theory. Weak lensing becomes large 
at small scales, and at large values of the perturbation wave- vector k perturbation series breaks down. The usual 
treatment of weak lensing therefore avoids cosmic perturbation theory by considering the small deflection angles 
of the photon trajectories. Another effect that has been estimated is i3-mode polarization from gravitational 
time delay (lo| and from sources proportional to second-order vector and tensor metric perturbations, which 
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are themselves generated from the product of scalar perturbations However, a full treatment of i?-mode 
polarization at second-order is still missing, even at the qualitative level, since previous second-order calculations 
considered the collision term [l^ and radiation transfer function 13] for unpolarized radiation only. In this 
paper, we derive the complete Boltzmann hierarchy at second order under the assumption that there are no 
first-order vector and tensor perturbations by extending the results of 12, 1^ to the photon polarization density 
matrix. This allows us to identify all sources of B-mode polarization at this order. The polarized equations are 
presented in a form suitable for numerical evaluation. Numerical results for B polarization will be presented in 



a follow-up paper IJ . Numerical results on non-gaussianity based on second-order Boltzmann equations have 



appeared recently in 1^ [l 

Most of the results of this paper have been obtained in the thesis work 17|. In the meantime, the polarized 
second-order Boltzmann equations have been derived independently in Our result is derived in a different 
formalism, allowing for an independent check of the results. We provide expressions for the Boltzmann hierarchy 
pertaining to the phase-space distribution functions not integrated over frequency, which have not been given 
explicitly before. In addition, we include a self-contained derivation of the polarized collision term from the 
quantum-mechanical time evolution of the photon density matrix, which differs from the collision term used 



18|. 



The outline of the paper is as follows. In the remainder of this section we set up our index conventions. Next, in 
Section In] we derive the Boltzmann equation from the quantum-mechanical time evolution of the photon density 
matrix. Its expansion to second order is presented in Sections IIIII and IIVI first for the propagation of polarized 
radiation (the "left-hand side" of the Boltzmann equation) , then for the collision term (the "right-hand side" ) . 
These sections are rather technical, and further technical details are collected in Appendix [XI The main result, 
the Boltzmann hierarchy for the second-order intensity and photon polarization phase-space distributions, is 
summarized in the separate Section [V] A full analysis of these equations is beyond the scope of this paper. 
However, in Section IVII we discuss the sources for _B-mode polarization, including a new source in the collision 
term that converts intensity directly into i?-mode polarization. We also analyze the tight-coupling regime. 
While we reproduce the presence of a second-order intensity quadrupole already found from the unpolarized 
equations [19|, we do not confirm the effect discussed in [20|, which is based on B-mode generation from an 
E-mode and intensity quadrupole in tight-coupling. Section [VIII summarizes our conclusions. 



A. Index and metric conventions 



General coordinate indices will be denoted by Greek letters n^v, . . . ranging from to 3, indices referring to 
tensors in the local inertial frame (tetrad frame) by capital Latin letters A,B,... = 0,1,2,3 from the beginning 
of the alphabet. Spatial indices, ranging from 1 to 3, in the tetrad frame are assigned Latin letters . . .. We 
also need small Latin letters a, 6, . . . = 1, 2 to denote the basis of polarization vectors. The signature convention 
for the space-time metric is (-1-,—,—,—). Spatial indices in the tetrad frame are contracted with the three- 
dimensional Euclidean metric and no distinction is made between upper and lower spatial indices. With this 
convention v^wa — v'^w'^ — v^w^ — v^w^ — v^Wi etc. In this paper we assume that the background universe 
is flat. We then also use Latin letters to denote the spatial general coordinate indices with the same 

convention regarding their contraction. In general, it will be clear from the context whether i, j, . . . refers to the 
tetrad or general coordinate system. Since confusion might arise for the momentum, we denote the covariant 
momentum dx^ /d\ by capital P^^, related to the momentum p"^ in the local inertial frame by = [e^]''p"^ 
with p^. 

The perturbed flat-space Robertson- Walker metric with conformal time denoted by rj and coordinates = 
(?7, a;* ) is parameterized as 

ds^ = {{1 + 2A) drj^ + 2B^ drjdx' - [(1 + 2D) + 2E,j] dx'dx' ) , (1) 

where 0(77) denotes the scale factor. The space-time dependent perturbations X — A, D, Bi, E^j will be expanded 

into first-order, second-order etc. terms according to X ~ X'^^^ + X'-^^ + We assume that the vector and 

tensor perturbations contained in B^, E^j are smaller than the scalar perturbations, so we formally treat them 
as second order. 
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We choose the conformal Newtonian gauge such that Bi is a transverse vector and Eij a transverse, traceless 
tensor. In this case, when there are only scalar perturbations at first order, we have = = 0. Then the 
tetrad components [e^]'* are given to second-order by 

[eo]° = ^ (l - A('^ - A(^) + ^A(i)' + ic/f ) ), 




HY = I {S., (l - D(^) - + ^D(i)^) - - (2) 

We denote by [e"^]^ the inverse of the tetrad, such that [e'^],^ [e^]'^ = Sg and [e'^]^[e^]'^ = S'^. The local inertia! 
frame can be parameterized in terms of the observer three- velocity Ui and a set of angles 6k, which define 
the orientation of the local inertial coordinate axes relative to those of x*. Above we have aligned the local 
coordinate axes with those of the general coordinate system to set the angles Ok = 0. In the following wc also 
choose the observer rest frame Ui = 0, which coincides to first order with another common frame choice Ui = Bi. 

II. BOLTZMANN EQUATION FOR THE POLARIZATION DENSITY MATRIX 

In this section we briefly review notation and definitions applying to photon polarization. We then derive an 

expression for the propagation and collision term in the Boltzmann equation for the polarization density matrix, 
which serves as the starting point for the expansion to second order in perturbations. 

A. Photon polarization phase-space distribution 

Wc assume that the polarized radiation ensemble can be described by a single-particle phase-space distribution 
matrix fij,u{x^, q^), such that e'^e"* fij,v{x^ , q^) <Pp/ (27r)^ denotes the number density of photons with momentum 
p and polarization e**. We regard f^^ as a function of the comoving momentum = ap'. The unperturbed 
Bose-Einstein distribution f^J is then independent of conformal time in the expanding homogeneous universe. 

The phase-space distribution is a Hermitian matrix, related to the expectation value {A^{x)Av{ii)) of the 
radiation field. We adopt Lorenz gauge A^.^ — for the photon field. It then follows that 

P^UA=o\q')=p''UA=o\q') = 0, (3) 
and that /^^ is parallel-transported in the absence of collisions. Thus 

^U. = c,Af]- (4) 

Here V/VX denotes the covariant derivative along a photon trajectory x^{X), and C';^y[/] is the collision term. 
The phase-space distribution distribution in the local inertial frame is related to f^i, by 

= [eAnenrf^'', (5) 

where [ba]^ are the space-time dependent tetrad vectors. 

The phase-space distribution matrix /^j^ is not unique, since Lorenz gauge allows the gauge transformations 
fnu ~^ ffiv + Oi^P^ -\- PuPfi with arbitrary a^, To obtain a physical distribution function, we decompose the 
photon four-momentum into 



P^' = E [eo]'' - [e,] V = K - n^), 



(6) 
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where = [cq]^ is the four-velocity of the locally inertial observer, E the energy of the photon as seen by this 
observer, and 

= (7) 

the photon three- momentum direction, which satisfies u^n^ = and n^n^ = — 1. We define 

Pi-iu = -ff^!^ + "m"!^ ~ nt^ny, (8) 
which projects on the components transverse to the observer velocity and photon direction: 

u'^'p^iv = u^p^i, = n^pp^ = n^p^i, = 0. (9) 
We now define the physical phase-space distribution matrix 

/m- =<V'iv', (10) 

which is orthogonal to P^, and n^, and contains no residual gauge ambiguity. The corresponding projected 
distribution function in the observer rest-frame (local inertial frame) is effectively a three-by-three matrix, since 
yoo _ yOi _ jio = Q jj^ frame. It is orthogonal to the photon propagation direction n* as expected for the 
tensor describing the two physical transverse polarizations. 
Applying V/VX to (HU]), we obtain from (g]) 

p/p/§^U'.' =p/p/c^'Afl (11) 

where we use that p^^, and /^^ are orthogonal to P^. This provides a closed set of equations for the physical 
phase-space distribution matrix provided the projected collision term C^!y[/] = p^^ pj^ C^i^i[f] on the right- 
hand side depends only on / rather than /. It will be convenient to choose a polarization basis consisting of 
two vectors e'^ (a = 1, 2) orthogonal to and , in terms of which 



r ^Y.f-^^'^'b- (12) 



B. Propagation of polarized photons 

In the following we shall consider the Boltzmann equation for fab in a polarization basis defined in the observer 
rest-frame. To this end, we multiply (fTTj) by e'^el" , insert (fT2|) . and calculate 

Before giving the result, we need to discuss the choice of polarization basis. If e** denotes the direction of the 



photon field-amplitude in Lorenz gauge, then is parallel-transported along the photon path [2l|. This also 
holds for the projected vector e'^ = p^ ^, i'^ in the sense that for the polarization vector corresponding to a fixed 
photon momentum p^^Ve^ /VX = 0. We could therefore use a basis of two parallel-transported polarization 
vectors e^. In this case, the corresponding vectors e;^ in the observer rest-frame, given by e{j' — [caI^Eq , depend 
on the space-time point x in addition to momentum q. The transverse polarization basis axes are continuously 
rotated as the photon propagates through space-time. 

It will be convenient to instead choose rigid basis vectors in the observer rest-frame that do not depend on 
space-time but only on momentum. In this case the covariant derivative along the path acting on a polarization 
vector is 



]^[e^]^,,P''e„Be^ + eaA^^, (14) 
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where the first term on the right-hand side arises, since is not parallel-transported in the rigid basis, and 
the second is due to the fact that the photon momentum also changes along the path. The conditions w^e^ — 
P^e^ = require = and n • = in the observer rest-frame. We adopt a spherical coordinate system with 

n — (sin 6 cos (p, sin 9 sin (p, cos 9) 
eg ~ (cos 9 cos (f, cos 9 sin (p, — sin 9) 

= (— siniy9,cosiy9,0) (15) 

the standard basis vectors on the sphere. The polarization basis is then taken to consist of the two circular 
polarization vectors 

e± ^ -^{eg ±ie^). (16) 



V2' 



We also identify — rj with conformal time. 

Evaluating (|13p results in the Boltzmann equation 



dU Idx^dU Idq^fdU, del* , * d4 r 
drj ^ P" dX dx, ^ d\\ dq^ ^ dq^ ^ Oq^ ^'^ 

+ [^kf-.u ^ {^a^efU + el,4fac) = C,, [/] , (17) 

with Cab[f] = ^a^l" ^ tf^[f] denoting the (projected) coUision term in the polarization basis. Its expression is 
given in ((39)) below. The terms to the left of the equality sign in the second line would vanish had we chosen 
a basis of parallel-transported polarization vectors. The extra terms in the rigid basis are equivalent to similar 



terms that appear in 22 1. 



For unpolarized radiation fab — Sab f, and the terms in ([T7| that depend on the polarization vectors explicitly 
vanish. This follows from 



since eak^l = ^ab, and from 



[euT,, (e«ef + ^eS) = Ca.e^'^ ([e% [ekf).^ = , (19) 



since [e*]^ [^kY — S^- Thus pT|) reduces to the standard equation for unpolarized radiation for diagonal phase- 
space density matrices, as should be the case. 

For the same reason, the explicitly polarization- vector dependent terms are at least of second order in per- 
turbations around the equilibrium distribution in the expanding homogeneous universe. This is due to the 
fact that dq^ /dX and [e*]^ are both first order in perturbations. Hence to first order we may set fab 

in the polarization- vector dependent terms equal to the unperturbed distribution f^^ ■ But the unperturbed 
distribution is diagonal, so the terms vanish (at first order) as shown above. 



C. Simplification at second order 

We now show that the terms to the left of the equality sign in the second line of (flT)) vanish even at second- 
order in perturbations, provided that there are no first-order vector and tensor perturbations. Thus, under 
these assumptions, there is no difference between the rigid and the parallel-transported polarization basis at 
second order. 

It follows from Section III Bl that the second-order contribution is the product of 

ea,e:Vi'^+44/ic' (20) 
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and 



(1) 



(0) 



B 



(1) 



dkB. 



(1) 



ni 



(5, 



(1) 



(i; 



(21) 



where the superscript in brackets indicates the order at which the expression is evaluated, and the elhpses denote 
terms proportional to qi or qk which vanish when contracted with the polarization vectors in (|20[) since q (x n. 
To obtain this expression we used the tetrads from (jAl[) . which do not assume a particular gauge choice. 



If there are no first-order vector and tensor perturbations, B^^^ and E}^' are zero in conformal Newtonian 
gauge, and expression (|2ip immediately vanishes, leading to the desired simplification. More generally, in an 
arbitrary gauge B^^'^ can be expressed as the gradient of a scalar function in the absence of vector modes, hence 
the curl of B^^^^ appearing in pT|) is zero. Likewise, E^j^ — {didj 
absence of first-order vector and tensor modes. Then, 



6ijd'^) E^^^ for some function E^-^'> in the 



which vanishes when contracted with the polarization vectors in (j20p . since ni£ai = 0. 



(22) 



D. Collision term 

To obtain the collision term Cab [/] we consider the quantum time evolution of the one-particle density matrix 
following the formalism developed in for neutrino flavour-mixing in a medium. The formalism was applied 
to photon polarization and Thomson scattering in 24] . A more general treatment elucidating some of the ap- 
proximations involved in the truncation of the hierarchy of n-particle density matrices implicit in this formalism 
can be found in [25|. 

In the local inertial frame with coordinates ^ the photon field operator is expanded in the form 

^(^) = E / {e-''-^'^aip)eaip) + e^^<aiip)e:ip)) . (23) 

We choose the two circular polarization vectors e± as basis vectors. The creation and annihilaton operators 
statisfy the standard commutation relation 

[aaip),alip')] = Sati27rf2p^6^'\p-p') ^ S^.Sip-p'). (24) 

The one-particle density matrix is defined by the expectation value {a\(j))aa{p')) ■ Spatial homogeneity implies 
that 

{al{p)aa{p')) = 5{p-p')pab{t,p). (25) 

We identify pabit,p) with the phase-space distribution function fi,a{x^,q^ = ap^). Indeed, since the number 
operator is 

^ = j:Jj^^-'aip)aa{p), (26) 

we obtain from (|^ 

N ={N)=V j d'piypit.p), (27) 

confirming the interpretation of pab{t,p) as phase-space polarization density matrix. The spatial dependence of 
fba{x^Tq^) can be neglected for the calculation of the collision term, since each scattering event is local on the 
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cosmological scales over which fba{x^,q^) varies. The flip in the order of polarization indices follows from the 
definitions (IT^ . and the fact that '^^aa{p)£a{p) is independent of the choice of polarization basis. 

The time evolution of the density matrix is obtained from the Heisenberg equation for the operator Dab{p) = 
al{p)aa{p)- Starting from 

J^Dab^i [H, Dab] , (28) 

going to the interaction picture and splitting the Hamiltonian into the free and interaction part Hj , we obtain 
to second order in the interaction [23| 

2^0(2^)3^(3) ^^pab{t,p) = i{[Hj{t),Dab{t,p)]} - £ dt' {[Hj{t - t'), [Hj{t), Dab{t,p)]]). (29) 

If Hi were the electron-photon interaction of quantum electrodynamics i?QED , we would have to expand to the 
fourth order in the interaction to recover the Compton scattering collision term. Instead we derive an effective 
Compton scattering interaction vertex assuming that the electron propagates freely between the two elementary 
electron-photon interactions in the Compton process. Thus we define Hi{t) through the relation 

i-i) J dtHiit) = J d^xd^vT{Hc^^^{x)H^^^{v)) (30) 

with the understanding that a pair of electron fields is contracted in the expression on the right-hand side. After 
a short calculation we obtain 24 1 

^ +p'" -q° -p") 



Hl{t) = J [dp][dp'][dq][dq'] {27Tf6'^^\q' + p' - q - p) , 

a,a' .s,s' 

X M{pa- qs p'a'; q's') a^, {q')al, {p')aa (p)Qs (?)■ (31) 

Here a, denote electron annihilation and creation operators, and [dp] = d'^p/ ((2Tr)^2p'^) is the phase-space 
integration measure. The matrix element for the 7(p,a) -I- e^{q,s) — > ■j{p',a') + e^{q',s') Compton scattering 
process reads 



C(pO /:y"^% ^a{p)+ Mp)/^4^ fa'iP') 

(q + P) —m^ [q — pY—mt 



u{q,s) (32) 



M{pa; qs -> p'a';q's') = u{q' , s') 
We note that 

M {pa; qs p'a'; q's') = M* {p'a'; q's' pa; qs). (33) 

To avoid confusion let us also note that in this subsection q stands for an electron momentum and not for the 
comoving photon momentum. 

The first-order term {[Hi{t), Dab{t,p)]) in involves the forward Compton scattering matrix element, 
and it is straightforward to show that this term vanishes. The second-order term is more complicated. It 
results in expectation values of four photon annihilation and creation operators, since the interaction generates 
correlations. To proceed we assume that n-particle correlations can be expressed in terms of one-particle 
correlations, such that, for example 

{al'{Q')aa{q)al,{p')ab{p)) = 5{q - p')Sab'{al,{q')ab{p)) + {al,{q')al,{p')aa{q)ab{p)) 

^ S{q-p')6ab'{al{q')ab{p)) + (4, (q')«a(g)) (4 (p')«fa(p)> + {al{q')ab{p)){al,{p')aa{q)) 

= 5{q ~ P')5{q' - p)pba'{p) [5ab' + Pab'{q)] + S{q - q')5{p - p')paa'{q)Pbb'{p)- (34) 



This amounts to the assumption that on average multi-particle correlations built up in a collision decay rapidly 
in the time interval before the next collision. The corresponding expressions for the electrons are simpler, since 
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we further assume that the electrons are unpolarized and that their phase-space density gdq) is sufficiently 
small for quadratic terms in g^, to be negligible. Thus 

(4(q')«.(9)4b')«r(p)> 5{q-p')5r>s{al,{q>r{p)) = 5{q ~ p')5{q' ~ p)5,r'5rs' \9e{q'). (35) 

Note that ge{q) is the density summed over both electron polarizations. After working out the expectation value 
of the second-order term in (|29p one ends up with the time integral 

t 

I '0 ^0 _0 



di'e±^*'(« +P ). (36) 







If the interaction time-scale is much shorter than the average time between collisions the upper limit may be 
taken to infinity and supplying the appropriate ie prescription, we obtain 

± ^ PV , „ + 7r5(g'° +p'°-q^-p\ (37) 

The imaginary principal-value term should be discarded, since it corresponds to a self-energy contribution. 
Putting everything together, substituting p^a — > fab in the last step, we obtain from ([^ 

2p°^/ab(p) = 2C„b[/], (38) 

where the collision term is given by 



= \ I (2,)4- (2.W° '^■')^''"''-' + " - - I" 



X <^9e{q)fx'u'{p) 5a\{S^b + Ub{p)) + 5^b{5aX + Ia\{p)) 

- 9e{q) Saxfuibip) + Sojbfaxip) {S X' u,' + f \' io' {p' ))^ ■ (39) 

Here we introduced the electron-spin averaged square of the Compton amplitude 



I M{p\- qs ^ p'A'; q's')M* [pcu; qs p'uj'; q's'). (40) 



The collision term (P^)) for the polarized phase-space density is the expression that must be used on the right- 
hand side of the Boltzmann equation ()17p . It takes an intuitive form with a gain and loss term and the expected 
Bose enhancement factors. Taking the trace in ab, and averaging the matrix clement over polarizations, we 
recover the standard unpolarized collision term. Eq. (pQ)) differs from [2j], where it is stated that the terms 
quadratic in the photon phase-space density cancel exactly in the evaluation of the double commutator in (j29|) . 



It also differs from the collision term used in [18|, which is based on |26i. The differences are located in the 
structure of the loss term from !26| and the Bose enhancement factors added in [l^ . The loss term in [26] is not 
derived as in the present paper but based on a certain ansatz, which is checked for initial and final pure photon 
polarization states, and then argued to hold in general due to the superposition principle. However, the loss 
term ansatz in [26| is non-linear in the phase-space distribution invalidating the superposition principle, and we 
suspect that this leads to the discrepancy with our result. Nevertheless, it turns out that the differences do not 
affect the final result in Section|V]below after the expansion to second order, at least for the frequency-integrated 



phase-space distributions considered in [18|. The reason for this is the simple polarization dependence of the 
Thomson scattering cross section and the fact that the terms quadratic in the photon phase-space densities will 
be seen to not contribute to the second-order equations for the frequency-integrated distributions. Differences 
between the present calculation and [18J from the form of the collision term would however be expected at the 
next order. 
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E. Fourier transformation and multipole expansion 



It is more convenient for the perturbation expansion to work with Fourier-transformed and multipole-expanded 
functions. We define 

(41) 

At second order we encounter products of functions, whose Fourier transform is a convohition. Below we use 
the short-hand notation 



(2^ y (2^(27r)=^'5(3Hfe-fci-fc2)A(fci)i?(fc2). 
For the multipole representation we write the comoving momentum OiS q — qn and then define 

oo / 



(42) 



.2^-1-1 

1=0 m=-l 



21 + 1 f 

fab,im{r],k,q) ^i^S^^— / dnYi%\{n) fab{r],k,qn). (44) 

Here Y(!^{ri) denotes the spin- weighted spherical harmonic. We collect the definitions and some basic relations 
for these functions in appendix lA 21 

We adopt the circular polarization basis (fT6|) such that under a rotation of the coordinate system around the 
direction of photon propagation with rotation angle A\l/ the polarization basis vectors transform according to 

<=±-e±^^*e,.±, (45) 

i.e. the circular polarization vectors e± have spin s = ±1 as they should. Since the polarization-basis independent 
phase-space distribution 



fV^J2C4fab (46) 



is invariant under basis rotations, it follows that f^^ and / arc spin-zero (s — 0) objects that do not 

transform, while 

/;^ = e±2^^*/±^. (47) 

Thus, has spin 2 and / ^ has spin —2. The corresponding values of s must be used in (|33]). 

Instead of the phase-space densities of the photon helicity states, one may also parameterize fab in terms of 
the four real Stokes parameters. The relation in the circular basis is 

The multipole decomposition for the Stokes parameter distribution functions fx reads 



fv,im = / dnvrMfvin), 



l2l + 1 f 

fEJm± lfB,lrn ^ \j J dQ F^^^* {n)±ifu{n)]. (49) 

The quantity // provides the photon density averaged over the two helicity states, and fv is related to the 
degree of circular polarization of the radiation plasma. We shall include fv in the set of second-order equations. 
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but since there are no sources of circular polarization in the standard cosmological scenario, it is usually of 
little interest. Our main concern are the ofF-diagonal components of the photon phase-space density, which are 
decomposed in (1491) into the E and B polarization modes. The conversion between the two sets of phase-space 
distributions follows from 



fx.lm = Ux:[ab]fabd 

Interpreting [ab] as a single index taking the values 
the matrix Ux;[ab] ^-J^d its inverse read 



, H — , — h m 



(50) 

this order, and with X — I,V, E, B, 



U 



X-[ab] 



0\ 








u. 



ab]:X 



V o-iiy 



/I 


-1 





o\ 


1 


1 














1 









1 





(51) 



We note the relations 



^[ab]-X — ^[bl]*X — 2t^X;[a6]- 



(52) 



In terms of multipoles the momentum derivative terms in the first line of the Boltzmann equation (|17p can 
be written in a simple form. First, from (jl6p we calculate 



T- 



^''dq' ' 9 tan 61 

where the upper (lower) sign holds for b = c = + (b = c = — ). Then, making use of (jA9|) . we obtain 
dfab , , de'^* J, , ^ _ dfgb ^ 2ie^, 



(53) 



^ak-g^fcb + ^bk'SZifo,'^ ~ ^ ■ + 



/+- 
gtanS ^ — / — ^ 



EH)' 

l.m 



47r dfab.lm 



21 + 1 



1 fab^lm / i 



V2 q 



(54) 



This form makes explicit that each term carries definite spin, such that s — for the diagonal elements and 
s = ±2 for the off-diagonals. The derivatives on the spin-weighted spherical harmonics can be easily taken using 
(IHOll . 



III. EXPANSION OF THE PHOTON PROPAGATION TERM TO SECOND ORDER 



We now turn to the expansion of the Boltzmann equation (|17p to second order in perturbations. Implementing 
the simplification of the polarization-dependent terms derived in Section III CI we obtain the first- and second- 
order equations 



d q^ d 
drj aE dx^ 



Jab 



1 dq' 



1' dfl 



(0) 
b 



q dq 



po 



Cabif] 



(1) 



(55) 



d q^ d 
drj aE dx^ 



(2) 
b 



pi 



1 dq* 



po 

(2) 



'-'Jab 1 


■ 1 dq^' 




dx* 







Cabif] 



(2) 



de: 



f-ak 



c Ai) 
Qqi Jcb 



de'y 



^bk 



a „i ac 



(56) 



Here we used that 
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at zeroth order in the perturbation expansion. In this section we keep the energy and momentum distinct so that 
the results also apply to the propagation of massive particles. For photons we may use E — \p\ — \q\/a — q/a and 
/q — to simplify the equations. Eq. (j55p reproduces the Boltzmann equation in the linear approximation 
with the familiar free-streaming term on the left-hand side. The Fourier transformation converts d /dx^ — ^ ifc' in 
the free-streaming terms. However, the second-order equation also contains products of two Fourier-transformed 
functions, which are to be interpreted as convolutions according to (021) • Thus, for instance. 



pi 



(1) C^fW 

'-'Jab 



pi 



(1) 



(2^ 



pt 



(1) 

{k-k')^k''fll\k'), 



(58) 



in the Fourier transform of (|56p . In this section we work out the multipole transformation of the left-hand side 
of (j56p . The more complicated transformation of the collision term is derived in Section ITVl 



A. Covariant momentum and momentum derivative 



The expression of the covariant momentum in terms of the comoving momentum required to evaluate (|56p is 
obtained from P'' = [cyij^p'^. Under the assumptions made in this paper (no first-order vector and tensor 
perturbations, conformal Newtonian gauge, observer frame, see Section [I A[) . we find 



a 



P' = l-D 



3^2 q'B, 
~2 



aE 



3D' 



n J^ki 



where the ellipses denote corrections of the third-order in perturbations. Hence, 



pi 



q^ 

aE 



(59) 
(60) 

(61) 



The change of comoving momentum dq^ /dX along the particle trajectory follows from the geodesic equation. 
We have 



dp^ ^ d{\e%P^ ^ p^p^ ^ (-Fi.VP^P") = [e\,.P^Pr 



dX 



Then 



dX 



1 dq' dq' da ^ a dp' ^ , P^P'' 



(62) 



(63) 



where He = a^^da/drj denotes the conformal Hubble parameter. The previous expression vanishes at zeroth 
order in the perturbations. Its perturbation expansion can be calculated from (1591) . (j60p and the explicit 
expressions for the inverse tetrad vectors. The first and second order terms required for (|56p read 



q 



J_d(f_ 
po rfA 

1 dq^ 
po dA 



(1) 



(2) 



aE 

q 

aE 



q'd'A^^^ -qi)^ 



(2) _ 9V ^(2) ^ oP ^(2) ^ ja^E^-q^) ^(2) 



q'B, 



qaE 



(64) 



(65) 



The dot denotes a derivative with respect to conformal time and 9' = d/dx\ The term proportional to He in 
the first line of (j65p vanishes for photons and massless propagating particles in general. 
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B. Multipole transformation and spherical basis 



The general procedure to obtain the multipole decomposition of ([55|) . ([SG]) is as follows. First we insert the 
representation (j43p for the phase-space distributions. Then the direction vector n and polarization vectors are 
written in terms of spherical harmonics according to 

Att Att 

= E^-V T^'"' '''''' " ^o'^^oo +E^2;.Vy>^2™, 

m m 

^+ = E/ ^™ \/~^^^™ ' ^ ~ ^™ Y''~3"-^1"i ' 

which defines (for m = 0,±1), xj/ = j5^^ and the trace-free tensors X2m (f*-"^ ~ 0, ±1,±2). Explicit 
expressions are provided in Appendix lA 31 The multiplication of these objects with Cartesian vectors and 
tensors, respectively, projects on the components of the corresponding vectors and tensors in the spherical 
basis. For vectors V and traceless symmetric tensors T we define the components in the spherical basis by 

3 

^[0] = ~2"^^^' "^'=^^1 ^ ±V2(ri3 T iT23) 

T[±2] = -^^1 - T22 T 2^Tl2). (67) 

Then 

X-imTi] = -ctmT[m] (no sum ovcr to) (68) 

with Qfo = |, a±i = and 02 = 1. At this point, we can use the product formula for the spin- weighted 
spherical harmonics (jA5[) to express any term in terms of a sum of single harmonics. The result of these 
manipulations is integrated with 



_4 121 + 1 



L^i'^j^^ j d^YiZin) (69) 

which projects (|55p . (j56p on the Im multipole component. The final step consists of transforming from the ah 
helicity polarization basis to the X — I^V,E,B components of the phase-space distribution matrix. 



C. Pree-streaming term 



We first consider the three space-time derivative terms in (|56p . which after Fourier transformation read 

1 (1) 



dt] ' aE ' 



po 



The multipole transformation of the time derivative is trivial since 

L 



drj 



(70) 



(71) 



For the transformation of the second term we follow the procedure described in Subsection IIII Bl The manip- 
ulations are the same as for the corresponding term in the first-order equation (|55p . and we discuss them here 
only to illustrate the general method. 



13 



Inserting the expansion of ~ qn^ and /^^^ in spherical harmonics gives 



iq k (2) _ 

aE ■'"^ aE 



47r, 

rn2\l 3 



li ,mi 



E 



in2 = — 1 
X 



47r 



2/i + l 



E E 

L=|(i-l| S,M=-L 



V2h 



hi L 

V2lTi V -s 



h 1 i 

TOi ni2 M 



(72) 



Applying the multipole transformation operator L from ([69|) to this expression sets L — I and M — m. 
Interchanging orders of summations according to 



oo ti L oo L + 1 L l\ 

EE E E=EE EE 

ii=0 mi=-li L=\l-i_-l\ M=-L L=0 li = \L-l\ M=-L mi=-h 



(73) 



yields the final result 
L 



iQ ■ k (2) 

Jab 



aE 



1 „7,[m2] '+1 'i 
nK ^ ^ 



/l 1 / 

-s -s 



h 1 I 

mi ni2 m 



fab,hmi(k)- 



(74) 



Here we used that the first Clebsch-Gordan coefficient sets S* = s to eliminate the sum over S. The second one 

implies m2 = m — mi. Recall that s takes the value when ab — H — h, and s — ±2 for ab — ±=p. Thus, in the 

first case only li = IzLl contribute to the sum, while for the off-diagonal terms li = I is also non-zero. Eq. (l74t 
reproduces the standard first-order free-streaming term, in which one usually aligns k with the three-direction 
implying fc^^^l — and k^^^ = ik, which simplifies the expression. 

The free-streaming term is diagonal in the circular polarization basis, but the equations for the two off- 
diagonal components are slightly different, which leads to a mixing of E and B polarization in the Stokes 
parameter basis. The difference arises from 



h 1 I 
2 2 



= (-ir 



h+i-i 



h 1 I 

-2 -2 



(75) 



i.e. when + 1 — Z is odd, which happens precisely for the terms with li — I present only for s — ±2. To express 
the equations in the IV EB basis in a compact form we introduce the matrices Hxy{1) with 



Hxy{1) = SxY (for I even). 



Hxvil) = 
















1 

















i 


Vo 





— i 


0/ 



(for I odd) 



(76) 



and define 

[-2 X = E,B 

Taking linear combinations of ([74| according to ([50l) we obtain in the IV EB basis 

1 , fm„l l+l h 



(77) 



L 



iq ■ fc .(2) 

aE 



E 



qk^' 



m2— — 1 
X 



aE 



E E 



li = \l-l\ mi=-li 



h 1 I 
Fx Fx 



h 1 / 

mi TO2 fn 



(78) 



Y 
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The sum over Y encodes the mixing between the E- and S-mode polarization. Since the H matrices are block- 
diagonal in IV and EB, and equal to the identity matrix in the IV sector, the sum is redundant for X — I,V. 
The equations are decoupled and identical ioT X = I and X — V. Nevertheless, the notation introduced above 
is convenient in order to present the results in the IV EB basis without having to resort multiple equations for 
the different cases. 

The third term in the list ([70)1 requires no further work, since using (pT|) 



(1) 



»g ■ fc2 .(1) 

aE 



(79) 



A convolution of the two mode momenta in the sense of (j42p is implied. The application of the L-operator 
gives as final result the expression ([71)) with fc^™^' — fc2'"^' and fj^\^jnii^) ~^ fab\imii^2), or the corresponding 
resuh dZS]) in the IVEB basis. 

The generation of B polarization from E polarization through free-streaming requires propagation through 
an inhomogeneous universe, and is thus a second-order effect, known as time-delay induced B polarization [id ]. 
The time-delay effect is contained in the above equations through the off-diagonal terms Hg^{±l) — —i. The 
relevant product of Clebsch-Gordan coefficients is 



1 I l \ I I 11 

2 2 /1 mi 1712 ITT- 



2m 



-'mi ,m— m2 



1(1 + 1) 



T\/2{l + lTm){l±m) 



TO2 = 
?7l2 = ±1 



(80) 



In first order in perturbation theory we can always align the mode vector k such that only TO2 — contributes. 
Then, using (1501) in ([75| . (1751) shows that EB mixing occurs only when to 7^ 0, which implies the well-known 
result that no B polarization is induced, when there are no vector or tensor perturbations. At second order ([79p 
contains a convolution over all wave-vectors, and the sum over TO2 always extends over TO2 = 0,±1. It follows 
from ([80[) that EB mixing occurs through free-streaming, when the first-order scalar perturbations A'^^^ or D*^^' 
do not vanish. 

To summarize the result of this subsection: the second-order space-time derivative terms (free-streaming 
terms) in the Boltzmann equation, Fourier- and multipole-transformed, are given in the IVEB basis by 



d_ 
dr] 



1 l + l h 

E E E ^ 

n2=-l = mi=-ii 



h 1 I 

Fx Fx 



h 1 I 

mi TO2 m 



Y^H*xy{ll + l-l) 



aE 



1 [7712] 



(81) 



D. Momentum-derivative terms 



We now turn to the multipole decomposition of the two terms involving dq^ /d\ in ([5^ . With the help of (IMl) 
and (|54l) the first one can be written as 



1 dq' 



1 (1) 



(1) 

b 



J rh ^ Hk , J ac 



de'l 



00 /l 

- (-1) E E H)'^ 

li—O mi— — /i 

aE 



Aw 



2/1 + 1 



— m-fciA(i)(fci)+i)(i)(fci) 
9 



— iklA^^^Uki) - -^ik\D^^\ki) 

q aE 



(82) 
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where we used that e± is orthogonal to q, and n - q = qn^ = q. Next we express n and the polarization vectors 
in terms of spherical harmonics according to (1661) to write 



in 



lm2 



m2 = — l 



(83) 



after taking the derivatives on the spin-weighted spherical harmonics using (jAlOl) in the second equation. The 
remaining steps are straightforward. We eliminate the products of spherical harmonics with (IA5I) and apply 
the multipole transformation operator ([M)) to obtain 



ihs of m 



l+l h 



ab 



-D''Hk,)q§-.fiZM+ E E E 



m2 = —l li=\l—l\ mi=— ii 



h 1 I 

-s -s 



aE 

q 



h 1 I 

mi 7712 m 



h 1 I 

(s - 1) -1 -s 



q ah 



(84) 



As for the free-streaming terms the equations for the off-diagonal terms are slightly different, which implies 
conversion of E into B polarization and vice versa. The last two lines in the previous equation, which originate 
from the derivative of the first-order photon perturbation with respect to the direction of the photon momentum, 
correspond precisely to the weak-lensing effect 0. If k{s) denotes the expression in curly brackets in the third 
line of the relation k{—2) = (— l)'i+i-' ^(2) holds, and because of the similarity with (1751) the same matrix 
HxY appears in the transformation to the Stokes parameters. The final result for this term in the IV EB basis 
reads 



L 



Ihs of m 



l+l h 



X 



-mki),§-/^\M+ E E E 



•i2 = -l /i=|i-l| mi=-ii 



h 1 I 

mi m2 m 



h 1 I 
Fx Fx 



aE\ 
1 ) 



h 1 I 
Fx-ll Fx 



X k' 



^A(')iki)-^D^'Hki) 
q aE 



)+^^'^^^^{fx\i-1Fx 

Y,H*xy{h + l-l)f^^l^S^2) 



(85) 



The other momentum-derivative term at second-order can be written as 



1 dq^ 



(2) J OJ-(O) 

q dq 



X + Y^^-Eg^^ 
q ■' q'^ 



(0) 



5abq 



dq 



(86) 



where X and Yi represent the q-independent and linear terms in in (|65l) . respectively. We also used that the 
unperturbed photon phase-space distribution is unpolarized. Since the only dependence on the direction of q 
in this term arises from the factors of g* in square brackets, it contributes only to Z = 0, 1, 2. In the quadratic 
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term we write 



p(2) 



5^(2) 



'^™2^[m2]-^2m2: 



(87) 



employing the definitions , and the tracelessness of Eij . The remainder of the calculation is straight- 
forward, resulting in 



L 



1 dq' 
aE 

q 



(2) 



1' dfit' 



q dq 



-D^^'>{k) + 2D^^'>{ki)D^^Hk2 



a„iE^^T^{k) Si2 '> Sabq 



(0) 



dq 



6ii 



(88) 



In the IV EB polarization basis the multipole transform of this term takes the same form with the replacement 
Sab Sxi- 

Our final result for the Boltzmann hierarchy for the multipole moments fx im^^) Stokes parameter 

phase-space densities at second order is given by the sum of (jHT]), ([85|l and ([88|) excluding the collision term 
that we consider in the following section. These expressions remain valid in the case of massive particles with 
mass M, for which 



E=\lAP + ^ 



(89) 



For photons the simplification aE /q = 1 can be applied and the term proportional to He in the second line of 
(l88ll vanishes. 



IV. EXPANSION OF THE COLLISION TERM TO SECOND ORDER 

In this section we compute the expansion of the collision term in the Boltzmann hierarchy for the multipole 
moments fximi'^)- This is done in two steps. First we expand (pQ]) to second order. Then we apply the 
operator (|69p that converts to equations for the multipole moments. Our treatment follows [l^I extended to the 
polarized phase-space distributions. 



A. Non-relativistic expansion 

The cosmic background photons that we see have mostly last scattered around the time of recombination, when 
the temperature of the universe was less than 1 eV. Polarization of the CMB is generated at this time or later. 
The electrons on which the photons scatter are therefore highly non-relativistic with thermal velocities 

(90) 

We therefore perform an expansion of the Compton scattering matrix element in the electron momentum and 
consider the expansion parameter (j90p of the same order as the cosmological perturbations. Note that in this 
subsection q and q' refer to the electron momentum and not the comoving photon momentum. 

The electrons are in local thermal equilibrium and sufficiently dilute to be described by the Maxwell- 
Boltzmann distribution 
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5pb 


(1) 

+ 


SXe 




. Pb . 









Here Te, Ve and rig denote the local electron temperature, bulk velocity, and number density of free electrons, 
i.e. electrons not bound in hydrogen or helium. If Xe denotes the ionization fraction and pi, the baryon density, 
then ng is given by 



(92) 



to first order in perturbations. A complete account of the collision term to second order therefore requires a 
calculation of the recombination history that goes beyond the homogeneous universe to obtain the perturbations 
in the ionization fraction. We refer to 27[ for a discussion of this issue. In our equations we keep rig as an 
overall factor without expanding it for the time being. 

The integral over q' in ((39)) is eliminated by the three- momentum delta- function, which sets q' ^ q + p p' ■ 
This allows us to expand 



7(9') ^ 9{p + Q-p') ^ g{q) 



[p - p'){q - rriev) (p-p'f , I f {p - p'){q - m^v) 



rripT 



2m.pT 



Sip" + -P'"- = S{p + E{q) -p'-Eip + q-p'))^S{p^ p') 



2 \ ?71eT 

{p-p')qd5{p~p') 
me dp' 



2mp 



dp' 



me 



[p - p'f d5{p - p') 1 ( {p-p')qV d'^5{p-p') 



dp 



/2 



where p = \p\t p' — \p'\- The expansion is based on the observation that p,p' ^ T while |q| 
that the difference of electron energies 



E{q)-E{q') 



•,'2 



q{p ~p') ip-p') 



,M2 



T, 



3/2 



2mp 2mp 



2mp 



^Te 



(93) 

(TOere)l/2 and 

(94) 



The terms neglected in (j93|) are therefore of third order in the expansion parameter (I90|) . Inserting these 
expansions into the collision term the zeroth-order terms cancel, so that the collision term begins at first order 
as it should be. It is therefore sufficient to expand the Compton matrix element to first order. The result of 
expanding ([5^ and (^0]) can be written in the form 



with 



^0. 



\M 



XX', 



XX' 



-SI. 



o 



£x{p) ■ e\,{p')el{p) ■e^,{p'), 



(95) 



(96) 



Si.xx'u^u;' ^ ex[P) ■ ex'iP ) i (P) + e^'(P ) ; 

p p 

<'{P')-P 



- elip) ■ e^. ip') { e\{p) lAl^^^ + e*,) [p') 



, ex{p) p' 



P 



p 



(97) 



and (Jt ~ 87rQ;^/(3TOg) the Thomson scattering cross section. At this point the integrand is polynomial in q 
except for ge{q) so that the integral over q in ([55]) can be expressed in terms of the moments of the electron 
distribution: 



dq 

(2^ 



ge{q) X {1; q'; q'q^] = x {l; meU*; meTpS^^ + mlvlvl} . 



(98) 
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B. Expansion of Ca.b[f] 



It is straightforward to insert the non-relativistic expansions discussed above into the colhsion term p9p and to 
perform the integrations over the incoming and scattered electron momentum. It is convenient to express the 
results in terms of the coefficient of the gain term in (|39|) and the difference of the gain and loss terms, given by 



(99) 



at zth order in the expansion. Note that while the difference of the gain and loss terms is linear in the phase- 
space distributions, the gain term contains quadratic terms. In the definition of G\y^^, we use the unexpanded 
distribution functions in the Bose enhancement factors. 

The expanded collision term can now be written in the form 



Cab[f] = ^neCTT I dp'p' I ^ 



-(1) -L ^(2) , „(2) , J2) , (2) , (2) 



(100) 



where denotes the solid angle of the scattered photon momentum vector p'. This expression includes the 
first-order term 



(1) _ o 
Cq(, — '^O.XX'uiui' 



'S{p - p') GL(^) + • (p - p') ^fc^ GL(«) 

op' 



(101) 



AA'i 



(summation over repeated photon polarization indices A, A',a;,a;' is understood), and the second-order term 
split into five contributions according to 



CAab = So,XX'uu- S{p - p') GL)^>,^^, 



(2) 



•'v,ab 



(2) 
'^Av.ab 



So,x...'vi^^.ip-p')''^^^rio) 



dp' 



'-'^XX'uiui' 



So,xx'..' vi'^ ■ {p - p') ^^-^ GL« + 5? 5{p - p') .(1)^ GL« 



dpf 



'XX'i 



(102) 
(103) 

(104) 



vv,ab 



So. 



1 r 



AA'. 



^ dH{p-p') rp^ 

Qp/2 ^^XX'uuj' 



+ "^l.AA'ww' IP P q^/ ^^XX'uuj' 



(105) 



(2) 
-K,ab 



[p-p'Y f dS{p-p') ^^^0) , dS{p-p') ^(0) d^5{p-^ ^0) 

0,XX'lulu' I — ^^\\'u„.>' ~ ^ ~ "-^AA',,,,,,' + TTT^ ^^XX'ujuj' 



y)2 

2me \ dp' 



+ S\^xX'u)U)' 



ip^p'Y 



-5{p-p')d^l^^,+T^ 



dS{p-p') ^^(0) 
dp' 



dp' 



(106) 



Due to the delta-functions the integral over p' can be performed after a few partial integrations. We also sum 
over polarizations and integrate over the solid angle, whenever possible. We define the integral operator 



2p 



1 fdn' 



(107) 



such that 



(108) 
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and work out the six terms separately. 
The first-order term yields 



(i)i 



1 f dn' 

An 



Sq. 



AA', 



\p=p 



(1) 



in-2n')GL\'l 



\p=p' 



dp' ^P=P 



(109) 



AA'aJoj' 



after partial integration. We note that 5*0 and SI depend on the direction of p' but not on its magnitude p' . 
The integral over the delta-function sets p' to pn' . The subscript "p = p'" means that fl]j{p') = fabipn') in 
the expressions (j99p for the gain and loss terms. The zeroth-order distribution function does not depend on the 
momentum direction and is unpolarized, hence 



\p=p 



= 0, 



AA'cjcj' 

Inserting the expression (IM|) for S'o,aa'( 



Pd^^^\P=P' 



(0) 



(110) 



into (|109p we next obtain 



einn')el:{n')fl^l{pn') - . (« _ «') el^{n')ei,in') p 



- e\ {n)e*J (n) 4? (n^^, (n') [SaxfH^ {p) + S^tf^'^ {p) 
This can be further simplified using 

y— .-(n')eUn') = y— 



(1) 



(0) 



dp 



kl Ik /li 

n n 



and e^(n)e"(n) = (^a^j to obtain the final result 



f^,l\p)^-Sabn-v^/>p 



dp 



dVl' 
47r 



2,. 



6;(n)6:^"(n) e*^,{n'y^,{n') f'^l,{pn') . (113) 



(111) 



(112) 



This expression is equivalent to the standard result for the first-order polarized collision term 28l. |29|. 

The second-order terms can be calculated in a similar way without further complications though the algebra 
gets lengthier, when the matrix element S\ is involved. We also note that the gain term alone, which contains 
quadratic terms in the phase-space distributions, appears only in c]^^;^ at second-order, but the simpler zeroth 
order expression 



gTx'..' ^^SaxS^bSx'., ff\p') [l + ff\p) 
is needed there. The result for the integrated second-order terms is 



P) 1 

'A,a6J 



I\c. 
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2 , df[°\p) 



,(2) 1 _ 



.ijkl 

f(i) 
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P ^ 1- Oab 
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4(n)e-(n)l.W^^)>^^ 



(118) 
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I\c 



47; 



m, 



dff\p) 
dp 



T 



dp"^ 



(119) 



Here S'™ jj,/ equals S\ xxi^ui' '^i^^ the polarization vectors stripped off, i.e. 



SI 



ijkl 



■jnil 

■' n 



(120) 



C. Fourier- and multipole transformation 



The final step in the derivation of the collision term consists in applying the multipole transformation operator 
(|69p . to perform the Fourier-transformation, and to convert the equations to the Stokes parameter basis. Taking 
into account the expansion of the prefactors 1 /P'^ in (fT7|) and the free electron density rig in (jlOOp , the right-hand 
side of the Boltzmann hierachy is of the form 
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at second order. In order to arrive at the last equality in each equation we used (|108p for the collision term, 
for the free electron density, and set E — p, since the collision term refers explicitly to photons. We 
also define the operation /[. . .] — L[I[. . .]]. The Fourier transformation of the collision term is trivial, since it 
does not contain spatial derivatives. Products of position-dependent functions simply turn into convolutions as 
indicated by the momentum argument. In the remainder of this subsection we neglect these arguments to avoid 
notational complications, but we restore them in the summary of Section [V] To complete the calculation of the 
Boltzmann hierarchy, it now remains to work out the multipole transformation of (I113p . and (jllSp - (jllQp . 

Inserting the multipole expansion p3)) for fyl^iipn') and applying the L-operator (|69|) the angular integrals 
can be expressed in terms of the matrices 
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where s = is implied for ah = H — h, and s = ±2 for ab = =p±. Since the polarization vectors are spin-1 

objects, the Q-matrices are non-zero only for I < 2. To transform to the Stokes parameter basis we use the 
matrices ([?T|) and define 



X,lm 



Ux,[ab]Q[ab],. 



(124) 



in analogy with ([50| for the phase-space distributions. In the IVEB basis the Q-matrices vanish for X = B 
for any value of /. The non- vanishing Q-matrices are given explicitly in (|A14p . The trace 
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is diagonal in the multipole indices and defines the w-symbols. Similarly in the IV EB basis 

tr {Qx,lrnQY,l'm') = Qx*iraQY,l'rn' = ^ ^XY ^U'^mm' (126) 

for 

4y = U*x-Aa^Ur.^.dAd- (127) 

Only very few of the w^V ^'"^ ^-^t zero. In particular, ujf^ — a;^^ = — -^/S/SO is the only off-diagonal term that 
couples to polarization. The other non- vanishing values are summarized in (jAlsp . 

After these preliminaries we turn to the explicit calculation beginning with the first-order term. With the 
definitions (I123P and (|125p it is straightforward to obtain from (|113p the expression 

- I l-fiZmiP) - Sa^Sn p'-^ + \ ^tU,f>Z.,aP)\ ■ (128) 



Note that the phase-space distributions fal\mip) helicity components of the electron bulk- velocity field 

''^e^im] depend on the Fourier mode vector k. This dependence is suppressed in this subsection as mentioned 

above. The transformation to the Stokes parameter basis requires the calculation of Ux,[ab] ^[^''^^af 1^°^ ^^^^ 
purpose we use 

2 UxM]Uial],x' 4'y'C^[A'c:'];y't/;,[A'.']/S™(^^) = 2 ^xx' u^%>5y,y f'^}^{p) 
= 2^^V/ff™(p) (129) 

and obtain 

i\c'^\ = \ {-/SL(rt - & v^l, p^I^ + J^yAUp)] (130) 

When inserted into ()12ip we reproduce the first-order collision term in the Boltzmann hierarchy for the polarized 
phase-space distributions in a notation similar to [sqI . The last term in brackets describes the generation of the 
i?-polarization quadrupole in Thomson scattering. 

(2) (2) 

The second-order terms can be calculated in a similar way. The first two, -^[c)^^^] and /[c^ ^j,], have the same 
structure as the first-order term and can be obtained from (|130p without additional work: 
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By far the most complicated expression to transform to multipole variables is the term /[c^^ ^j,] in (jll7p . In 
terms of the Q-matrices and w-coefficients introduced before we find 
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We have made the sums over Y, Z expUcit here. While providing a closed expression for any X = I,V, E, B, this 
result is not very transparent. Recalling that the Q-matrices and w-coefhcients vanish for l,li, L > 2 and noting 
that the Clebsch-Gordan coefficients are non-zero only if the angular momenta differ by no more than one, we see 
that I[c^^y x] vanishes when ^ > 3. For any particular X the sums can be worked out explicitly at the expense 
of introducing explicit values of the Clebsch-Gordan coefficients. We give the corresponding simpler expressions 
in our summary of the Boltzmann hierachy in Section |Vl The last two terms to be converted to the multipole 

(2) (2) 

representation, I[CyJ ^^ij] and /[c^^j], are relatively simple, since they depend only on the unperturbed phase- 

(2) 

space densities. In particular, ^[c]^^^] does not contain any perturbation variables, hence angular dependence, 
since it arises from the non-relativistic expansion. The result is 
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Here and above in (jl33l) we expressed the result directly in the Stokes parameter basis. The result in the circular 
polarization basis is obtained by omitting the sums over Y,Z and replacing H*pq{. . .) — > 5pQ (any P,Q); by 
substituting 6xi Sab, X ^ ab, Y ^ X'uj' , Fx —s, Fy — >■ —s' (with s chosen according to the value of ab 



and s' according to A'w'), as weh as WxV 

mm - firm . 
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V. BOLTZMANN HIERACHY AT SECOND ORDER - SUMMARY OF EQUATIONS 

At this point we can return to our original notation and express the phase-space densities in terms of the 

comoving momentum q = ap, using pdfx{p)/dp = qdfx/dq. In the following wc leave away the photon 
momentum argument q on the phase-space densities but restore the Fourier mode momentum writing fx,im{ki). 
As mentioned above by taking the four values of X separately, we can evaluate the angular momentum sums 
over I, I' etc., and obtain a more explicit form of the Boltzmann hierarchy. Wc; summarize the second-order 
equations in this section. For convenience we recall the first-order equations in the absence of first-order vector 
and tensor modes in the present notation: 

= I'^l { - /ll(fe) + ^ioA%{k) - Sn1^vi]l^ik) + 5i2^ - ^/6/gL(fe)) } (136) 



d_ 
dt] 
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{ - /gL(fc) - (/i'L(fc) - v^/gL(fc)) } (138) 



y/bU^) + ^{Ti)f^B\i^,)^Sk)k^"^^^D^l^ + z/£L,(fc)fc[-]i?5kU = l«l { - /SL(fe)} (139) 
Here we introduced the abbreviation 

k = -n'^^^ara < (140) 

for the collision rate. Furthermore, here and below a summation over m2 = 0,±1 is implicitly understood in 
terms containing the index m2, and mi is equal to m — m2. We also introduce the coupling coefficients 
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as well as 



(142) 



Note that we may choose k such that it points into the three-direction, in which case = and the first-order 
equations become particularly simple. 

We now present our main result, the Boltzmann hierachy for the second-order perturbations to the polarized 
phase-space densities. Recall that the equations are given in conformal Newtonian gauge for a comoving and 
aligned observer (Ui = 0,6i = 0) under the assumptions of vanishing first-order vector and tensor modes 
{bI^^ = E^f =0). In the equations given below wc keep terms involving the first-order perturbations /jj^]„j 
of the i?-polarization density to display their structure. Of course, under the above assumptions there is no 
B polarization in first order, so fg\^ vanishes, and the corresponding terms can be neglected in numerical 
evaluations. The equations read: 
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(146) 



These are the dynamical equations for the second-order photon variables. The source terms depend on products 
of first-order perturbations as well as on the second-order perturbations A^'^\ D^'^\ -^[m]' -^\m\ ^'^ t'^^ metric 



and to the bulk electron velocity j^^j . To close the system of equations, these quantities must be determined 
from the second-order Einstein and fluid equations. 



At this point it seems appropriate to compare our results to those given in [18|. We already mentioned 
that the collision term in 18| takes a different form before expansion of the phase-space distributions around 
the equilibrium distributions, but that these structural differences drop out at second order, at least for the 
frequency-integrated equations. The derivation of the expanded equations in fisf follows a different method 
from the one empolyed in the present paper by first considering the collision term in the electron rest frame, and 
then performing the boost to the frame, in which the electron fluid moves with bulk velocity v^- In contrast, we 
work directly in this frame adopting the Maxwell-Boltzmann distribution (j9ip for the electrons. Both methods 
should give the same results, since the Lorentz non-covariance of the shifted Maxwell-Boltzmann distribution 
is a higher-order effect. For a detailed comparison we note that only the frequency-integrated equations for the 



quanitites A^j^(77, fc) defined in (|148p below are given explicitly in [l8| and that the contribution from ck in 



(I135P is neglected. The integrated equations can be obtained from the above by applying the substitution rules 
(|149l) . After doing this we find that the structure of the equations is in complete agreement but we observe 
differences in the following terms: the octupole collision source term for i?-mode polarization (the (5(3 term in 
our (|145p ) has different numerical coefficients (this is corrected in the arXiv version of [l^); in the B-mode 
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equation (our (|146l) ) the coupling coefficient ^^A™ differs from our corresponding D^^^^ ^ in the collision term 
and second-order Liouville operator, and the terms corresponding to to the last line before the equality sign in 



(I146P are missing 31 1. 



VI. DISCUSSION 

While a numerical or even qualitative evaluation of the second-order Boltzmann hierarchy is beyond the scope 
of the present paper, we briefly discuss the sources of i?-polarization contained in the equations, and the tight- 
coupling limit. Before proceeding to the discussion of the collision term we note the different I dependences in 



the weak-lensing and gravitational time-delay terms 32|, which we identify as the product terms of A^^\ D'^^ 



a mode momentum ki or k2 , and fx\m '^^ the left-hand side of the Boltzmann equations (|143p - (|146|) . While 
lensing of X = I,V,E,B on itself is proportional to I for large I, since R^'^K^'^ cx I (for large I), the 
corresponding time-delay effect is only of order 1, since C*'',!?^'' oc 1. In contrast, for conversion of E- into 
B-polarization and vice versa, weak lensing and time delay are effects of the same order, and both coefficients 
involved, if^^' and are only of order l/l for large I. 

A. S-mode polarization from scattering 

There are two sources of -B-mode polarization in the photon propagation terms on the left-hand side of the 
Boltzmann equations. A well-known mechanism is the generation of B polarization when polarized radiation 
propagates through an inhomogeneous universe, usually referred to as the weak-lensing effect. It appears first 
at second-order and is contained in the terms involving the product of the metric perturbation A^^"^ or D^^^ 
with the first-order i?-mode distribution f^\^ in the last two lines before the equality sign in (|146p . 

B-mode polarization is further generated in the presence of vector or tensor metric perturbations. Around pho- 
ton decoupling Thomson scattering generates the vector and tensor components of the i?-polarization quadrupole 
which is subsequently partially converted to _B-polarization through free-streaming. In the present scenario we 
assume that there are no first-order vector or tensor metric perturbations. In the absence of any primordial vec- 
tor or tensor perturbations, they will still be generated at second-order, however. B-mode polarization induced 



by these second-order perturbations through free-streaming has been estimated in [11|. The effect turns out 
to be relatively small, though comparable to the weak-lesning effect in the small Z-region of the BB anisotropy 
spectrum. 

The full second-order Boltzmann equations exhibit further sources for B polarization through the collision 
term, which are absent in the first-order equation (jl39p . which contains only the damping term —fslmi^) 
the right-hand side. The second-order collision term in (|146|) contains products of the electron velocity and 
first-order intensity and i?-mode perturbations. Of particular interest is the term 

^<L(^^^ (4 + 2) fi:LH^^)D"J,rn. (147) 

which can generate a i?-mode quadrupole directly from the intensity quadrupole rather than indirectly through 
iJ-polarization. A numerical analysis of the B-polarization generated from this term will be presented in [141] . 

B. Tight-coupling limit 

We now examine the second-order equations in the regime where the electrons and photons are strongly coupled 
by Thomson scattering. For the following discussion, we are not interested in the frequency dependence of the 
photon distribution functions and integrate over q. We define the frequency-integrated multipoles 
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In the fluid description of photon radiation Aj"gQ equals the fractional perturbations of the photon number 

density, and Aj'^j'^ = Av^^^-^ is related to the bulk velocity of the photon fluid. 

Using partial integration, derivatives on photon distributions can be eliminated, resulting in the following 
substitution rules in the Boltzmann equations in Section IVl 
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JX,lm ~^ ^X,lm 

1^ ^ -4A(,"L. (149) 
The only term to which these rules cannot be applied is the cx contribution from (jl35p to the collision term 

(2) 

for /) QQ, which contains non-linear terms in the photon distribution. Inserting the Bose-Einstein distribution 
for the zeroth-order /j"-' to calculate (|148p for this term, we find 



me me 



oq rrif. oq'^ rrie 



But in the strongly coupled electron-photon plasma the electron and photon temperatures coincide, so this term 
makes no contribution to the frequency-integrated Boltzmann equations. 

In the tight-coupling regime the collision rate |k| is larger than any other scale of interest. The collision term 
drives the system to equilibrium, which makes the left-hand sides of the Boltzmann equations small. Thus 
the Boltzmann equations can be satisfied only, if the coefficients of |k| in the collision term on the right-hand 
side nearly vanish. At leading order in the expansion in 1/|k| this enforces a number of relations among the 
perturbation variables. 

Looking at the first-order equations in Section |V] we immediately find ^y\„-^ = ~ ^ '^'-"^ ^'^^ 

intensity equation ()136|) has no collision term for / = 0, so the intensity monopolc is unconstrained in the 
tight-coupling limit. For the dipole ? = I, we obtain the familiar relation 

Ai!L(fc)=4^;W„j(fe), (151) 

which implies that the bulk velocities of the photon and electron plasma are equal. (The continuity equation 
for the electron fiuid yields the same relation.) Continuing with the quadrupoles, we obtain from p36p . (|138p 
the equations 



2 a(i) - -^A^') fl52^ 

which imply ^^i\m ~ 2m ~ Likewise, all higher multipoles vanish. It follows that there is no polarization 
in the tight-coupling limit, as expected, and only the intensity monopole and dipole are unsuppressed. 

We now consider the second-order equations in the tight-coupling regime. It is straightforward to see that 

f2l (21 

as in first order, circular and B polarization vanish, AyJ^^^ — A^ = 0, as well as the multipoles higher than 
I ~ 2 for / and E. The collision term for the intensity monopole is no longer zero, but vanishes at leading 
order in the tight-coupling expansion after inserting the relation (|15ip . so the monopole is again unconstrained. 
Setting I — 1 in (|143p we find the tight-coupling relation 

AS^(fc) = 4 (v^^l^{k) + .(;)„,(fci)A«o(fc2)) , (153) 
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which is similar to psip but contains a term quadratic in the first-order perturbations. Finally, we examine 
the quadrupoles. For / = 2, we can write (|143p . (|145p in the form 



with 
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Setting the right-hand sides of (|154p to zero yields the tight-coupling relations 

M^L(^) = 10We,[mi](fcl)We,[m2](fc2)C,;;f„, 



(155) 



(156) 



We therefore find that there is no polarization in tight-coupling at second order. However, contrary to the first 
order, there exists a non-vanishing intensity quadrupolc quadratic in the electron bulk velocity, as expected. In 
cartesian components and before Fourier transformation, Eq. (jl56p corresponds to 



(157) 



in agreement with 19[, where this result has been obtained from the unpolarized Boltzmann hierarchy. (A 



factor of two difference arises due to the different convention for expanding quantities X to second order.) It 
follows from the above that the size of the quadrupole is not modified when the full polarized set of equations 
is employed, since the i?-polarization quadrupole vanishes in tight coupling. 

Our results are at variance, however, with (2o| . where it has been found that the tight-coupling intensity 
quadrupole provides a large source for S-mode polarization. The argument is based on an incomplete expression 
for the i?-polarization source term (|155l) . Since the authors of [io^ did not have the Boltzmann equations for E- 
and S-polarization at second order available, the source term without the product of first-order perturbations 
was used, 
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The source term was further simplified using the second of the relations (I152p for the second-order modes as done 
after the second arrow above, which implies the assumption that there exists ^^-polarization in tight-coupling 
in contradiction with (|156p . We conclude that in this case it is clearly important that the full second-order 
polarized equations are used. Then it follows from (|155p that P'^^ = 0, and thus there is no source term in 
tight coupling that would yield E- and therefore _B-polarization from the line-of-sight solutions of (|154[) . The 
large effect reported in 20] is therefore absent. Polarization is only generated, also at second-order, once the 
scattering rate drops sufficiently so that corrections to tight coupling become relevant. 



VII. CONCLUSION 



In this paper we derived the complete Boltzmann hierarchy for the polarized photon phase-space distributions 
at second order in conformal Newtonian gauge and in the local observer rest frame under the assumption that 
vector and tensor perturbations are formally of second-order. This assumption is well-motivated by the fact 
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that our primary aim is to study the B-mode polarization and non-gaussianity induced at second-order, when 
the primordial sources are small. A first analysis shows that the B-mode collision term contains new sources 
that involve the intensity of the perturbation rather its i?-polarization. In tight-coupling we obtain the intensity 
quadrupole found earlier from the unpolarized Boltzmann hierarchy but no i?-mode polarization. The equations 
presented here set the stage for their numerical evaluation, which we plan to present in a subsequent paper. 

Acknowledgements. We thank P. Fischer for collaboration during an early stage of this work, C. Pitrou for 



correspondence on Ref. [18| and helpful suggestions, and K. Klingmiiller and Y.Y.Y. Wong for comments on the 
manuscript. M.B. wishes to thank the theory groups at the University of Zurich and CERN for their hospitality 
during 2008 and 2009, when part of this work was performed. This work is supported in part by the Gottfried 
Wilhelm Leibniz programme of the Deutsche Forschungsgemeinschaft. 

Appendix A: Summary of definitions 
1. Tetrad components 

The tetrad components before specifiying confornial Newtonian gauge are given to second order by 
a 

-Ie^e^ (A1) 

Quantities without superscript are expanded according io X ~ X'^^ + X^^^ + . . .. The one simplification that 
has been made is that we set to zero the angles Ok , which defines the orientation of the local inertial coordinate 
axes relative to those of a;*. The expressions in conformal Newtonian gauge adopted in this paper are given 
in ©. 



2. Spin-weighted spherical harmonics 

The spin- weighted spherical harmonics are defined for ? > |s| and \m\ < I by 



+ (-1)'"''""+™ e*™'^ cot2'^+"-" (A2) 



such that for s = the standard spherical harmonics are recovered. carries spin s, since under a rotation 
of the coordinate system with angle IS.^ is transforms as 

r'L = e-^*r4. (A3) 

For any given s the spin- weighted spherical harmonics define a complete set of functions on the sphere obeying 
the orthogonality relations 

' dnYiZYi^,,,=5u.5,r,rn'. (A4) 
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Under complex conjugation Yi^* = (— l)™^''y/_^- A product of two spin-weighted spherical harmonics can be 
combined to a single one using 



hrn^ Z^^ y/4TT{2l + 1) \m1m2mj \ -Si -S2 S J ^ ^ 



The summation ranges are restricted by the triangular equation for the Clebsch-Gordan coefficients 

( M ^0 if 1^2-^11 <Z<^i+^2, (A6) 

Y mi 7712 m J 

and 

M = if 771^^7711+ m2. (A7) 
TTJl TO2 ml 



This implies in particular s — si + S2 in (|A5p . Furthermore, the Clebsch-Gordan coefficients satisfy the following 
relation: 



= (-1)'!+'^-' " ] . (A8) 

mi TO2 ml \ —mi — m2 — m ' 



The spin-raising and -lowering operators are defined by 
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scot 9, 



We then have 



where 
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d i d „ , , „, 

-7^ + — -scot 0. A9 
o9 sm dtp 



If = v/('Ts)(;±.s + l). (All) 

3. Unit vector in the spherical basis 

The coefficients and x^m defined in ()66p which express n* and n' n-' in terms of spherical harmonics are given 
explicitly by: 

eo = I oj ^±i = _L 1^^, I , (A12) 

. /-I 0\ / Tl\ ^ / 1 0^ 

X20 = -10 X2,±i = * X2.±2 = -= -1 I . (A13) 
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4. Q matrices and uj coefficients 



The non-vanishing Q-matrices introduced in (|124p read in the IV EB basis: 
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The non-zero traces (|126p are: 
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